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Abstract 

Background: Interferon-gamma (IFN-gamma) and interleukin-6 (IL-6) are multifunctional cytokines that regulate 
immune responses, cell proliferation, and tumour development and progression, which frequently have functionally 
opposing roles. The cellular responses to both cytokines are activated via the Janus kinase/signal transducer and 
activator of transcription (JAK/STAT) pathway. During the past 10 years, the crosstalk mechanism between the 
IFN-gamma and IL-6 pathways has been studied widely and several biological hypotheses have been proposed, but 
the kinetics and detailed crosstalk mechanism remain unclear. 

Results: Using established mathematical models and new experimental observations of the crosstalk between the 
IFN-gamma and IL-6 pathways, we constructed a new crosstalk model that considers three possible crosstalk levels: 
(1) the competition between STATl and STATS for common receptor docking sites; (2) the mutual negative 
regulation between SOCSl and S0CS3; and (3) the negative regulatory effects of the formation of STATl/3 
heterodimers. A number of simulations were tested to explore the consequences of cross-regulation between the 
two pathways. The simulation results agreed well with the experimental data, thereby demonstrating the 
effectiveness and correctness of the model. 

Conclusion: In this study, we developed a crosstalk model of the IFN-gamma and IL-6 pathways to theoretically 
investigate their cross-regulation mechanism. The simulation experiments showed the importance of the three 
crosstalk levels between the two pathways. In particular, the unbalanced competition between STATl and STAT3 
for IFNR and gpl30 led to preferential activation of IFN-gamma and IL-6, while at the same time the formation of 
STATl/3 heterodimers enhanced preferential signal transduction by sequestering a fraction of the activated STATs. 
The model provided a good explanation of the experimental observations and provided insights that may inform 
further research to facilitate a better understanding of the cross-regulation mechanism between the two pathways. 



Background 

The Janus kinase/ signal transducer and activator of tran- 
scription (JAK/STAT) pathway can be activated by a 
wide range of cytokines. The binding of cytokines to 
their receptors drives receptor dimerization and phos- 
phorylation, which leads to the recruitment, activation 
and dimerization of STATs [1-5]. STAT homodimers are 
translocated to the nucleus, where they regulate the ex- 
pression of target genes [6,7] . IFN-gamma was first iden- 
tified in a mitogen-activated lymphocyte supernatant 
based on its distinctive antiviral activity. It is widely 
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involved in protection against tumour development and 
cancer immunoediting [4]. IFN-gamma triggers the pro- 
longed activation of the transcription factor STATl via 
the IFN-gamma receptor (IFNR) and JAK, which 
induces target gene expression by binding to the 
gamma-activated sequences (GAS) in the promoters of 
these genes [6]. Thus, the activation of the IFN-gamma/ 
JAK/STATl pathway can prevent the expansion of many 
different normal and neoplastic cells [8]. However, IL-6 
also has important roles in triggering the acute phase re- 
sponse of the body to injury or inflammation [2]. The re- 
ceptor involved with the recognition of IL-6 can be 
subdivided into the non-signalling part, i.e., glycoprotein 
80 (gp80), and the signalling part, i.e., glycoprotein 130 
(gpl30). IL-6 attaches to gp80 and drives IL6-gp80 
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complex binding to the gpl30-JAK complex, which forms 
the IL6-gp80-gpl30-JAK complex. gpl30 associates with 
JAK and it is tyrosine-phosphorylated in response to IL-6 
stimulation, which leads to the activation of the JAK/ 
STATS pathway and MAPK (mitogen-activated protein 
kinase) cascades [2,9]. IL-6 activates target genes involved 
with differentiation, survival, apoptosis and proliferation, 
and it plays important roles in pro- and anti-inflammatory 
functions, acute phase and immune responses in the or- 
ganism, and tumour progression [2,6]. Three types of 
negative regulators are involved in the regulation of the 
IFN-gamma and IL-6 pathways: the suppressor of cyto- 
kine signalling (SOCS), SH2 domain-containing tyrosine 
phosphates 2 (SHP-2) and various cytoplasmic and nu- 
clear protein tyrosine phosphates (PPs) [10,11]. SOCSl 
and SOCS3 are induced by JAK/STAT pathway and they 
bind to the activated receptors of IFN-gamma and IL-6, 
respectively, which negatively regulate the signal transduc- 
tion of IFN-gamma and IL-6 [12,13]. SHP-2 acts as a 
phosphatase in activated receptor complexes of IFN- 
gamma and IL-6 and it negatively regulates the activation 
of STATs [14,15]. PPl dephosphorylates STAT*s (* repre- 
sents species-activation states in this study) in the cyto- 
plasm, while PP2 dephosphorylates STAT*s in the 
nucleus, which results in STATs being returned to the 
cytosol, thereby influence the activation of STATs [16,17]. 

Previous studies have indicated that IFN-gamma and IL- 
6 have opposing roles in cell proliferation, apoptotic death 
and inflammation, which are closely related to the specific 
patterns and duration of STAT activation after their stimu- 
lation [6]. IFN-gamma mainly phosphorylates STATl, 
which has many pro-inflammatory effects [8,18,19]. By 
contrast, IL-6 is a potent activator of STAT3, which contri- 
butes to its anti-inflammatory functions [2,20,21]. How- 
ever, the detailed molecular mechanism leading to the 
unbalanced activation of STATs after IFN-gamma and IL-6 
stimulation remains unclear. Qing et al. suggested that 
tyrosine 419 in the IFN-gamma receptor subunit 1 
(IFNGRl) is required for the activation of both STATl and 
STAT3 [19]. In response to IL-6, STAT3 binds to phos- 
phorylated (p)YXXQ motifs (Y767,Y814,Y905 and Y915) of 
gpl30, whereas STATl is recruited to a more restricted 
consensus sequence pYXPQ (Y905 and Y915) in gpl30 
[2]. In addition, researchers have provided some interesting 
experimental results using STAT-deficient cells. Qing et al. 
showed that the activation of STAT3 in response to IFN- 
gamma was much stronger and more prolonged in 
STATl-null cells than wild-type cells [19]. Costa-Pereira 
et al. showed that IL-6 mediated an IFN-gamma-like re- 
sponse in mouse embryo fibroblasts lacking STAT3, in- 
cluding the prolonged activation of STATl and it 
promoted the induction of multiple IFN-gamma-inducible 
genes [20]. However, Regis et al. reported that the activa- 
tion of STATl in human neoplastic T lymphocytes after 



IFN-gamma stimulation was generally unaffected by 
STAT3 silencing [22]. Ho et al. argued that STAT3 did not 
affect the activation of STATl [23]. Bluyssen et al. reported 
that pre-treatment of EC with IFN-gamma significantly 
reduced the activation of STAT3 after induction by IL-6, 
but without affecting the total amounts of STAT3 [18]. 
However, Kaur et al. reported that the activation of STATl 
by IFN-gamma was mainly unaffected after pre-treatment 
with IL-6 or other gpl30-related cytokines in SH-SY5Y 
human neuroblastoma cells [24]. These experimental 
results indicate that STATl and STAT3 may have common 
binding sites within IFN-gamma and IL-6 receptors, while 
the activation of STAT3 may depend on the concentration 
of STATl and vice versa. Moreover, the interactions be- 
tween IFN-gamma and IL-6 signals are not symmetric. 

Thyrell et al. also reported that IFN-alpha could influ- 
ence the signal response of IL-6 in multiple myeloma, 
which resulted in a decrease in STAT3 homodimer DNA- 
binding activity and a shift from STAT3 homodimers to 
STATl/3 heterodimers [25]. Herrero et al. observed that 
pre-treatment with IFN-gamma could affect the signal re- 
sponse of IL-10 in macrophages, which caused the IL-10 
mediated STAT activation pattern to switch from STAT3 
homodimers to STATl/3 heterodimers [26]. Therefore, 
changes between STAT homodimers and STATl/3 het- 
erodimers may represent a biologically relevant approach 
to determining the crosstalk between IFN-gamma and IL- 
6 pathways. However, how the formation of STATl/3 het- 
erodimers regulates the interactions between IFN-gamma 
and IL-6 signals is not fully understood. 

Systems biology modelling generally aims to find fairly 
plausible mechanistic models that include all of the relevant 
key processes in a biochemical system. It is considered to be 
a powerful analytical approach to understanding the essen- 
tial mechanisms of the physiological functions of normal tis- 
sues and pathological progression during complex diseases 
[27]. The first systems biology model of the IFN-gamma/ 
JAK/STATl pathway was developed by Yamada et al. [28]. 
They modelled activation of the JAK/STATl pathway in re- 
sponse to IFN-gamma and analysed the effects of the feed- 
back loop regulated by SOCSl. Later, Zi et al. conducted a 
multi-parametric sensitivity analysis of the model produced 
by Yamada et al. and indicated that the concentrations of 
SOCSl, nuclear phosphatase PP2 and cytoplasmic STATl, 
as well as some of the reaction steps that affected those con- 
centrations, were the most sensitive to perturbation [29]. 
The first model of the IL-6/JAK/STAT3 pathway was pro- 
duced by Singh et al. [30]. Recentiy, Moya et al. proposed an 
updated model of IL-6 and IL-10 signalling via JAK/STAT 
and ERK-C/EBPp activation [31]. The model was used to in- 
vestigate dynamical features of the system such as the activ- 
ity ratio of JAK/STAT and ERK-C/EBPp with different 
stimulation levels of IL-6 and IL-10. The dynamic beha- 
viours of some individual molecules, such as STATs and 
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SOCSs, in the IFN-gamma and IL-6 pathways were investi- 
gated in previous studies [28,30,31], but signalling of the 
crosstalk during signal transduction by IFN-gamma and IL- 
6 has still not been modelled. 

In this study, we developed a crosstalk model of the 
IFN-gamma and IL-6 pathways by combining previously 
established mathematical models and by comprehen- 
sively analyzing the interactions between the two path- 
ways. The model considered three possible levels of 
crosstalk between the two pathways: (1) the competition 
between STATl and STAT3 for IFNR and gpl30; (2) the 
mutual negative regulation between IFN-gamma and IL- 
6 via the regulators SOCSl and SOCS3; and (3) the 
restrictive effects of the formation of STATl/3 heterodi- 
mers on the activation of the transcription factors 
STATl and STAT3. We considered a number of proto- 
cols where cells were stimulated by IFN-gamma and/or 
IL-6. The simulation results showed that the model pro- 
vided a good explanation of the experimental observa- 
tions and it provided new insights that could inform 
further research to facilitate a better understanding of 
the cross-regulation between the IFN-gamma and IL-6 
pathways. 

Results 

Model description 

Based on the model of the IFN-gamma/ JAK/STATl 
pathway produced by Yamada et al. [28] and the model 
of the IL-6/10/JAK/STAT3 pathway produced by Moya 
et al. [31], we established a crosstalk model of the IFN- 
gamma and IL-6 pathways. A schematic diagram of the 
model is shown in Additional file 1: Figure SI. In this 
model, the components of the two previous mathemat- 
ical models, their structures and most of the parameters 
were left unchanged. For simplicity, we specified that 
SHP-2 could repress the activated receptors of IFN- 
gamma and IL-6, while PPl and PP2 (PPX and PPN in 
Yamada et al.'s work) could dephosphorylate STATl and 
STATS in the cytoplasm and the nucleus, respectively. 
We removed any reactions and components that were 
not connected with IFN-gamma and IL-6 signalling, 
such as IL-10. Sixteen new reactions were added based 
on the possible mechanisms of cross-regulation between 
IFN-gamma and IL-6. 

The structure of the STATl and STATS proteins contains 
an oligomerization domain, a coiled-coil domain, a DNA- 
binding domain, a linker domain, an SH2 domain and a 
transactivation domain [32,33]. The recruitment of STATl 
and STAT3 to the activated receptor complexes is known 
to be mediated by their SH2 domains and phosphorylation 
of the receptor tyrosine motifs is required [34-36]. After re- 
ceptor binding, the STATs are phosphorylated on a single 
tyrosine residue (Tyr 701 in STATl and Tyr705 in STAT3) 
[37,38]. Many experimental observations have shown that 



STATl and STAT3 may combine with the same docking 
sites in IFNR and gpl30 [2,19]. Thus, we hypothesized that 
STATl and STAT3 might compete for the same phos- 
phorylated docking sites in IFNR and gpl30, via their SH2 
domains in our model. After STATl and STAT3 combine 
with the activated receptors complexes via IFN-gamma 
and/or IL-6, they are phosphorylated and disassociate from 
the receptors. Based on these considerations, the new bio- 
chemical reactions (Nl) - (N3) and (N4) - (N6) (see "New 
biochemical reactions added to the crosstalk model" in 
Additional file 1) were added to our model to simulate the 
activation of STAT3 after IFN-gamma stimulation and the 
activation of STATl after IL-6 stimulation, respectively. 

The kinetic parameters of these new reactions are im- 
portant because they reflect the properties of the bio- 
logical system. Wiederkehr-Adam et al. indicated that 
the SH2 domain of STATl had a much higher affinity 
for the phosphotyrosine 419 motif in IFNGRl than that 
in STAT3 [39]. After IL-6 stimulation, STAT3 binds to 
the flexible pYXXQ motifs in gpl30, whereas STATl is 
recruited to the more restricted consensus sequence of 
pYXPQ in gpl30 [2]. Based on these observations, we 
hypothesized the unbalanced competitive binding of 
STATl and STAT3 with IFNR and gpl30 after IFN- 
gamma and IL-6 stimulation, respectively. Additional file 1: 
Tables SI -S3 show that the main effector of IL-6 sig- 
nalling, STAT3, had a higher affinity for gpl30 than 
STATl. Similarly, the main effector of IFN-gamma sig- 
nalling, STATl, had a higher affinity for IFNR than 
STAT3. 

SHP-2 and SOCS combine to regulate signal transduc- 
tion by IFN-gamma and IL-6 [1,12,40]. SOCSl inhibits 
the JAK/STAT pathway by binding to the activation loop 
of JAK via its SH2 domain [12]. SOCS3 can also bind to 
JAK [41]. SOCSl and SHP-2 combine with different sites 
in the receptor complexes of IFN-gamma [28] . However, 
SOCS3 and SHP-2 may have similar binding specificities. 
Experiments have suggested that SOCS3 and SHP2 may 
compete for same site (tyrosine Y759) in gpl30 after IL- 
6 stimulation [42,43]. In our model, SOCSl and SHP-2 
were capable of binding to the receptor complex of IFN- 
gamma without mutual interference, whereas SOCS3 
and SHP-2 could competitively bind to the receptor 
complex. Distinct genes belonging to the SOCS family 
are induced as immediate early genes (lEGs) down- 
stream of different STATs and they can inhibit STAT ac- 
tivation in a classical negative feedback loop [44]. It is 
generally recognized that SOCSl has an important role 
in modulating IFN-gamma signalling [45], whereas 
SOCS3 mainly affects IL-6 signalling [46]. It is also well 
known that the STATl and STAT3 homodimers are dir- 
ect transcription factors of the JAK/STAT pathway, 
which play important roles in signal transduction during 
IFN and gpl30 receptor signalling [6]. However, 
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Bluyssen et al. reported that SOCS3 could also be 
induced after IFN-gamma stimulation in EC and that it 
could inhibit signal transduction by IL-6 [18]. Qing et al. 
observed that the STATl and STAT3 homodimers could 
both be induced after IFN-gamma stimulation in MEFs, 
which both bound to the same GAS element in the 
SOCS3 promoter [19]. The sequence of this GAS elem- 
ent is conserved in mice, rats and humans [47]. It was 
shown that STATS activation was much stronger and 
more prolonged in STATl-null cells, and that SOCS3 
was strongly induced in wild-type and STATl-null cells, 
while the levels of SOCS3 mRNA were greatly increased 
in STAT3-null cells [19]. Thus, it is speculated that 
STATl homodimers might also promote the transcrip- 
tion of SOCS3 in the same way as STAT3 homodimers. 
However, no experimental evidence indicates that 
STAT3 homodimers can combine with the promoter re- 
gion of SOCSl. Thus, our model does not regard STAT3 
homodimers as an efficient transcription factor for 
SOCSl. We added equation (N7) to our model to simu- 
late the transcription of SOCS3 mRNA after its induc- 
tion by STATl homodimers in the nucleus, which is 
represented as (STAT1N*)2. 

Thyrell et al. reported that IFN-alpha could affect the sig- 
nal response of IL-6 in multiple myeloma, which resulted 
in a decrease in STAT3 homodimer DNA-binding activity 
and a shift from STATS homodimers to STATl/3 heterodi- 
mers [25]. Herrero et al. showed that pre-treatment with 
IFN-gamma could affect the signal response of IL-10 in 
macrophages, causing the IL-lO-mediated activation pat- 
tern to switch from STATS homodimers to STATl/3 het- 
erodimers [26]. These experimental results showed that 
STATl/S heterodimers play important roles in the cross- 
talk between different cytokines. The activation of STATs 
after cytokine stimulation led to the formation of STAT 
homo- and heterodimers [2,6]. 

Haan et al. reported that IL-6 stimulation of primary 
human macrophages led to a different distribution of 
STAT dimer species in the cytosol and nucleus. In particu- 
lar, they showed that STATl/S heterodimers were present 
in the cytosol and nucleus [48]. The size of STATs exceeds 
90 kDa, which is far beyond the exclusion limit of the nu- 
clear pores, so STATs need to be translocated actively into 
the nucleus [2,6]. Tyrosine phosphorylation is not neces- 
sarily required for STAT nuclear translocation [49]. Basic 
residues (Lys410 and Arg41S) are known to contribute to 
the nuclear localization signals (NLSs) of dimeric STATl 
[50]. PPl dephosphorylates STAT*s in the cytoplasm, while 
PP2 dephosphorylates STAT*s in the nucleus, which leads 
to STATs being returned to the cytosol [16,17,51]. It was 
postulated that the nuclear export signals (NESs) in the 
DNA-binding domain of STATl is comprised of the resi- 
dues 399-410 [52]. The formation and dephosphorylation 
of STATl and STATS homodimers in the cytosol and 



nucleus were modelled by Yamada et al. and Moya et al., 
respectively. Their models also described the translocation 
of STATl and STATS homodimers from the cytosol to the 
nucleus, and the export of dephosphorylated STAT mono- 
mer from the nucleus to the cytoplasm [28,S1]. In our 
model, biochemical reactions (N8) - (N9) were added to 
simulate the formation of STATl/3 heterodimers in the 
cytoplasm and nuclei. For simplicity, it was supposed that 
only STATl* and STATS* could form STATl/S heterodi- 
mers. Biochemical reaction (NIO) was added to simulate 
the STATl/3 heterodimer translocation process from the 
cytoplasm to the nucleus based on the translocation of 
STAT homodimers. It was also assumed that STATl/S het- 
erodimers could be dephosphorylated by PPl and PP2, 
which resulted in STAT export from the nucleus to the 
cytoplasm. Biochemical reactions (N11)-(N16) were added 
to simulate this process. The specific biological role of 
STATl/S heterodimers remains obscure [6], so we did not 
regard the STATl/3 heterodimers as efficient transcription 
factors in our model. 

To construct the crosstalk model, we merged the com- 
mon components from previous models, such as SHP-2, 
PPl and PP2. JAKl and JAK2 are two species in the JAK 
family, which play important roles in the signal responses 
of IFN-gamma and IL-6 [40]. In previous studies, JAKl 
and JAK2 were treated as JAK for simplicity [28,31]. JAKl 
and JAK2 can combine with the receptors of IFN-gamma 
and IL-6, but the signal transduction activities of IFN- 
gamma and IL-6 may relate to specific types of JAKs. The 
internal membrane proximal regions of JAKl and JAK2, 
which responds to the IFN-gamma signal, bind the IFN- 
gamma receptor subunits IFNGRl and IFNGR2 [53]. Dur- 
ing IL-6 signalling, JAKl and JAK2 are activated via the 
conserved membrane proximal-binding domain of the 
receptors, and JAKl plays a major role in the signal re- 
sponse to IL-6 [54]. Moreover, Guschin et al. suggested 
that, although JAK2 was activated, it could not mediate 
the efficient activation of STATl and STATS after IL-6 
stimulation in the absence of JAKl [55]. Thus, JAK2 may 
share redundant functions with JAKl. In our model, we 
used two different JAK species for IFN-gamma and IL-6 
receptors, respectively. 

Our model contains two main components: IL-6 
signalling via the JAK/STATS pathway and IFN-gamma 
signalling via the JAK/STATl pathway. There are multi- 
level interactions between the two pathways. The model 
contains 108 species, 192 kinetic parameters and 119 
reactions, of which 103 reactions are based on previous 
models whereas 16 reactions were new. To make the ex- 
perimental results easier to compare, if not specified 
otherwise, the concentrations of STATl* and STATS* 
were the sum of the concentrations of all species con- 
taining activated STATl and STATS, respectively, in- 
cluding their monomers and dimers. 
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Responses of the crosstalk model to separate IFN-gamma 
and IL-6 stimulation 

First, we stimulated the model with IFN-gamma (0.1 nM) 
for 12 h and found that STATl* reached its maximum 
concentration (350 nM) within about 1 h, before it 
decreased rapidly due to the feedback inhibition of SOCSl 
and SHP2. It finally reached a new steady state (60 nM) 
after about 6 h. STATS was activated to reach its max- 
imum concentration (60 nM) within about Ih and it 
decreased rapidly to the control level after 2 h, whereas 
the activation of STATl was much stronger than STATS 
after IFN-gamma stimulation (see Figure lA). The signal 
transduction profiles of these molecules were consistent 
with previous experiment results, although there were 
some differences in the signal strength and duration [18]. 



Next, we stimulated the model with IL-6 (0.1 nM) for 
12 h and found that STATS* reached its maximum concen- 
tration (S50 nM) within about 0.5 h, before it decreased 
rapidly due to feedback inhibition from SOCSS and SHP2. 
It reached a new steady state (40 nM) after about 6 h. 
STATl was activated and reached its maximum concen- 
tration (100 nM) within about 0.5 h, before it decreased 
rapidly to the control level after 1.5 h, whereas the activa- 
tion of STATS was much stronger than STATl after IL-6 
stimulation (Figure IB). The signal transduction profiles 
of these molecules agreed with experiment results [20]. 
Next, same kinetic affinities were set for IFNR and gplSO 
in STATl and STATS, respectively. As a result, IFN- 
gamma and IL-6 stimulation caused similar strong activa- 
tion of STATl, STATS, SOCSl and SOCSS (Figure IC-D). 



STAT3* 

STATl* 

(STAT3N»)2 

STAT1N«-STAT3N» 

— (STAT1N»)2 
SOSC3 

SOSCl 

SHP2 




5 6 7 
Time(h) 

Figure 1 Simulated time courses of species In the crosstalk model with continuous exposure to IFN-gamma (0.1 nM) or IL-6 (0.1 nM). 

(A) IFN-gamma stimulation induced much stronger activation of STATl than STATS. (B) IL-6 simulation induced much stronger activation of 
STATS than STATl. (C) In a balanced competition model, IFN-gamma stimulation induced similar activation of STATl and STATS. (D) In a balanced 
competition model, IL-5 stimulation induced similar activation of STATl and STATS. (E) With no formation of STAT1/S heterodimers, IFN-gamma 
stimulation induced more STATl and STATS homodimer formation than normal conditions. (F) With no formation of STATl/S heterodimers, IL-6 
stimulation induced higher STATl and STATS homodimer formation than normal conditions. In the Figure, SHP2 represents the monomers of 
inactive SHP2, which did not combine with other molecules or complexes. 
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The balanced activation of STATl and STATS after IFN- 
gamma and IL-6 stimulation did not agree with previous 
experimental observations. These results demonstrated 
the validity of our unbalanced competition model. 

We also investigated the signal transduction profiles of 
STAT homo- and heterodimers in the nucleus after IFN- 
gamma and IL-6 stimulation, separately. After continuous 
stimulation with IFN-gamma (0.1 nM) for 12 h, (STATIN*) 
2 reached its maximum concentration (65 nM) within 
about 1 h and it maintained a new steady state (5nM) after 
6 h, whereas (STAT3N*)2 only reached its maximal con- 
centration (1 nM) after 1 h (Figure lA). By contrast, IL-6 
stimulation (0.1 nM) for 12 h made the (STAT3N')2 level 
reach its maximum concentration (90 nM) within about 
0.5 h and it reached a new steady state (5 nM) after about 
6 h, whereas (STAT1N*)2 only reached its maximal concen- 
tration (1 nM) after 0.5 h (Figure IB). Our results con- 
firmed the experimental observations of Haan et al. who 
showed that IL-6 stimulation led to STATS homodimers 
predominating in the nucleus [48] . These results suggested 
that IFN-gamma and IL-6 signalling preferentially activate 
nuclear STAT homodimers. For the STATl/S heterodimers 
(STATIN'-STATSN*) in the nucleus, however, both IFN- 
gamma and IL-6 could induce a similar concentration/ 
strength, which reached its maximum concentration (50 
nM) in about 0.5-1 h (Figure 2A-B). IFN-gamma and IL-6 
could both activate STATl and STATS, but fewer STATl* 
and STATS* molecules were sequestered by STATl/S het- 
erodimers, so its transcriptional activation function was 
repressed. When the formation of STATl/S heterodimers 
was blocked, the maximum concentrations of (STAT1N)*2 
and (STATSN)*2 both increased to about 15 nM after IFN- 
gamma or IL-6 stimulation (Figure lE-F). Thus, the forma- 
tion of STATl/S heterodimers enhanced the preferential 
signal transduction of IFN-gamma and IL-6. 

SHP-2 knockout simulations were also performed to 
characterize the effects of SHP-2. As shown in Figure 2A-B, 
knocking out SHP-2 enhanced the signal responses of 
IFN-gamma and IL-6, which agreed with previous experi- 
mental observations. You et al. showed that IFN-gamma 
could induce a higher signal response in SHP-2 null cells 
[14]. Schapter et al. also reported that over-expression of 
an inactive SHP-2 mutant in HepG2 cells enhanced STAT 
activation after IL-6 stimulation [15]. After IFN-gamma or 
IL-6 stimulation, however, the JAK/STAT pathway exhib- 
ited different features to those when knocking out SHP-2. 
Without SHP-2, IFN-gamma stimulation induced higher 
levels of STATl* and STATS* (about 150%) than that in 
normal conditions (Figure 2A). By contrast, IL-6 stimula- 
tion induced rapid increases in STATl* and STATS*, the 
maximum concentrations of which reached 900 nM and 
500 nM, respectively, which was about three times higher 
than that in normal conditions. After IL-6 stimulation, we 
also observed that SOCSS reached a peak value about 75 



nM at 2 h, which inhibited signal transduction by IL-6 
and quickly caused the concentrations of STATl* and 
STATS* to drop to normal levels after S h (Figure 2B). 
Knockout simulations were also performed for SOCSl 
and SOCSS. As shown in Figure 2C-D, knocking out 
SOCSl enhanced the activation of STATl after IFN- 
gamma stimulation, while knocking out SOCSS enhanced 
the activation of STATS after IL-6 stimulation. Our simu- 
lation results agreed with previous experimental observa- 
tions. Brysha et al. demonstrated that in vitro and in vivo 
hepatocytes lacking SOCS-1 exhibited prolonged activa- 
tion of STATl after IFN-gamma stimulation, which corre- 
lated with the dramatically increased sensitivity to the 
toxic effects of IFN-gamma [IS]. Niwa et al. reported that 
inhibition of SOCSS expression enhanced the activation of 
STATS and cell growth [56]. After IFN-gamma or IL-6 
stimulation, however, the JAK/STAT pathway exhibited 
different features when knocking out SOCSl or SOCSS. 
Without SOCSl, IFN-gamma stimulation induced higher 
levels of STATl* and STATS* (about 250%) compared 
with those in normal conditions (Figure 2C). Without 
SOCSS, however, IL-6 stimulation induced increases in 
STATl* and STATS*, the maximum concentrations of 
which reached 700 nM and SOO nM, respectively, which 
were about double those in normal conditions. After IL-6 
stimulation, we also observed that SHP-2 dropped to a 
low level of about 80 nM at 1 h, which attenuated signal 
transduction by IL-6 and caused the concentrations of 
STATl* and STATS* to fall slowly after S h (Figure 2D). 
Thus, SOCSl and SOCSS had dffferent interaction pat- 
terns with SHP-2. SOCSl and SHP-2 synergistically regu- 
lated signal transduction by IFN-gamma. Knocking out 
SOCSl or SHP-2 enhanced the integral activation of 
STATl induced by IFN-gamma stimulation. By contrast, 
SOCSS and SHP-2 regulated signal transduction by IL-6 
in a more complementary manner. Knocking out SHP-2 
alone enhanced the rapid response of the IL-6 signal, due 
to a compensatory increase in SOCSS. Knocking out 
SOCSS also led to lower levels of SHP-2, which caused a 
slow decline in STATl* and STATS* S h after IL-6 stimu- 
lation. Simulations of the combined knockout of SHP-2 
and SOCS were performed to characterize their joint 
effects on IFN-gamma and IL-6 stimulations. First, we sti- 
mulated the SHP-2 and SOCSl combined knockout 
model with IFN-gamma (0.1 nM) for 12 h and found that 
STATl* reached its maximum concentration (960 nM) 
within about 2 h while STATS* reached its maximum con- 
centration (900 nM) within about S h (Figure 2E). IFN- 
gamma stimulation induced a similar strength of STATl* 
and STATS* in SHP-2 and SOCSl combined knockout 
conditions. We then stimulated the SHP-2 and SOCSS 
combined knockout model using IL-6 (0.1 nM) for 12 h 
and found that STATl* and STATS* both quickly reached 
their maximal concentration of 950 nM and 980 nM, 
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Figure 2 Simulated time courses of species in the crosstalk model with continuous exposure to IFN-gamma (0.1 nM) or IL-6 (0.1 nM) 
after knocking out SHP2 and/or SOCSs. (A) Without SHP-2, IFN-gamma stimulation induced higher levels of STATl* and STAT3* (about 150%) 
than normal conditions. (B) Without SHP-2, IL-5 stimulation induced rapid increases in STATl* and STATS*. (C) Without SOCSl, IFN-gamma 
stimulation induced higher levels of STATl* and STAT3* (about 250%) than normal conditions. (D) Without SOCSS, IL-6 stimulation induced rapid 
increases in STATl* and STATS* the maximum concentrations of which reached 900 nM and 500 nM, respectively. (E) With SHP-2 and SOCSl 
combined knockout, IFN-gamma stimulation produced a similar strength of STATl* and STATS*, which both rapidly reached their maximum 
concentrations of 960 nM and 900 nM, respectively. (F) With SHP-2 and SOCSS combined knockout, IL-6 produced similar strengths of STATl* 
and STATS* which both rapidly reached their maximum concentrations of 950 nM and 980 nM, respectively. (G) With SHP-2 and SOCSl 
combined knockout, the level of (IFN-R-JAK)2* increased significantly after IFN-gamma stimulation. (H) With SHP-2 and SOCSS combined 
knockout, the level of (IL6-gp80-gplS0-JAK)2* increased significantly after IL-6 stimulation. 
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respectively, within about 1 h (Figure 2F). IL-6 stimulation 
also induced similar strengths of STATl* and STATS* in 
SHP-2 and SOCS3 combined knockout conditions. Thus, 
the combined knockout of SHP-2 and SOCSs abolished 
the preferential activation of IFN-gamma and IL-6. The 
unbalance competition between STATl and STATS was 
not related directly to SHP-2 and SOCSs, but SHP-2 and 
SOCSs combined with the activated receptor complexes 
and inhibited signal transduction via the JAK/STAT path- 
way [2,6,12]. Therefore, we deduced that SHP-2 and 
SOCSs could limit the concentration of active receptor 
complexes, which indirectly affected the preferential acti- 
vation of IFN-gamma and IL-6. Therefore, we investigated 
the signal transduction profiles of the activated receptor 
complexes in response to IFN-gamma and IL-6 with and 
without knocking out SHP-2 and/or SOCSs. Without any 
knockout, (IFN-R-JAK)2* reached its maximum concen- 
tration (0.1 nM) in about 0.5 h after IFN-gamma, before 
decreasing rapidly. After knocking out SHP-2, the level of 
(IFN-R-JAK)2* increased by about 120% compared with 
that in normal conditions. Without SOCSl, (IFN-R-JAK) 
2* increased rapidly and reached a new steady state (0.S5 
nM) after 2 h, whereas the combined knockout of SHP-2 
and SOCSl caused the level of (IFN-R-JAK)2* to increase 
significantly, reaching 4.5 nM in 12 h, which was about 
forty time as high as that in normal conditions (Figure 2G). 
Without any knockout and with IL-6, (IL6-gp80-gpl30- 
JAK)2* reached its maximum concentration (0.15 nM) 
within about 0.25 h, before decreasing rapidly. After 
knocking out SHP-2, (IL6-gp80-gpl30-JAK)2* increased 
rapidly and it reached its maximum concentration of 1.4 
nM in 0.25 h, which was about nine times that in normal 
conditions, although it quickly returned to a normal level 
after 0.5 h. With SOCSS knock out, the (IL6-gp80-gpl30- 
JAK)2* level increased and reached a new steady state (0.15 
nM) after 1 h. With the combined knockout of SHP-2 and 
SOCSS, the levels of (IL6-gp80-gplS0-JAK)2* increased 
significantly and reached a new steady state (5 nM) after 
1 h, which was about 35 times that in normal conditions 
(Figure 2H). The simulation results demonstrated that with 
the SHP-2 and SOCSs combined knockout, the levels of 
(IFN-R-JAK)2* and (IL6-gp80-gplS0-JAK)2* increased sig- 
nificandy after IFN-gamma and IL-6 stimulation. STATl 
and STATS competed for the same motifs in IFNR and 
gplSO, but there was sufficient (IFN-R-JAK)2* and (IL6- 
gp80-gplS0-JAK)2*, so the preferential activations of IFN- 
gamma and IL-6 were abolished. These simulated observa- 
tions still await further experimental verification. 

Responses of the crosstalk model after disrupting STAT1 
and STATS 

The effect of STATS on signal transduction via the JAK/ 
STAT pathway was analyzed by varying the initial concen- 
tration of STATS in a range of 0-2000 nM. We found that 



changing the STATS level did not significantly affect the 
state of STATl* after IFN-gamma stimulation (Figure SA), 
which was consistent with previous experimental observa- 
tions [22]. By contrast, the level of STATl* was clearly 
affected by the initial STATS concentration in response to 
IL-6. In particular, when STATS was knocked out, STATl 
was more phosphorylated and for longer, so STATl* 
reached its maximum concentration (180 nM) in about 
1 h, which was about double that in normal conditions. Fi- 
nally, it reached a new steady state (50 nM) after about 
7 h (Figure SB). This was consistent with previous experi- 
mental results, although there were some differences in 
the signal strength and duration [20]. The different signal 
responses to IFN-gamma and IL-6 during STATS disrup- 
tion may explain why IL-6, but not IFN-gamma, could 
trigger apoptosis and inhibit the in vivo growth of human 
malignant T cells after knocking out STATS [22]. Next, we 
analyzed the effect of STATl on signal transduction via 
the JAK/STAT pathway by varying the initial concentra- 
tion of STATl in a range of 0-2000nM. We found that 
changing the initial concentration of STATl did not sig- 
nificantly affect the level of STATS* after IL-6 stimulation 
(Figure SC). By contrast, the level of STATl dramatically 
affected the status of STATS* after IFN-gamma stimula- 
tion. When we knocked out STATl in our model, IFN- 
gamma stimulation also led to much stronger activation 
of STATS, which caused a significant increase in the levels 
of STATS*. It finally reached a new steady state (220 nM) 
after 1 h, which was about three times that in normal 
conditions (Figure SD). Our simulation results were con- 
sistent with previous experimental observations [19]. 
Therefore, IFN-gamma and IL-6 signalling could mutually 
switch in the conditions of STATl or STATS knockout, 
which agreed very well with previous experimental obser- 
vations [19,20]. 

After IFN-gamma and IL-6 stimulation, the recruit- 
ment of STATl and STATS to the activated receptor 
complexes directly affected their phosphorylation, which 
has important roles in signal transduction by IFN- 
gamma and IL-6 [S4-S6]. In our unbalanced competition 
model, STATl and STATS had different affinities for 
IFNR and gplSO, so we deduced that disrupting STATl 
and STATS may have different effects on the recruit- 
ment of STATl and STATS. Next, we investigated the 
effect of changing the initial concentration of STATl 
and STATS on the associations of STATs with activated 
receptor complexes in response to IFN-gamma and IL-6. 
Our simulation results demonstrated that changing the 
concentration of STATS had little effect on the forma- 
tion of (IFN-R-JAK)2*-STAT1 (Figure 3E), whereas alter- 
ing the STATl level significantly affected the formation 
of (IFN-R-JAK)2*-STATS (Figure 3F) after IFN-gamma 
stimulation. Our simulations also showed that the for- 
mation of (IL6-gp80-gplS0-JAK)2*-STATS was almost 
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Figure 3 Dependency of JAK/STAT pathway activation on the initial concentrations of STATs. (A) The activation of ST ATI was mainly 
independent of tiie concentration of STAT3 after IFN-gamma simulation. (B) The concentration of STAT3 significantly affected the activation of 
ST ATI after IL-6 simulation. (C) The activation of STAT3 was mainly independent of the concentration of STAT1 after IL-6 simulation. (D) The 
concentration of ST ATI significantly affected the activation of STAT3 after IFN-gamma simulation. (E) Recruitment from STAT1 to (IFN-R-JAK)2* 
was mainly independent of the concentration of STAT3 after IFN-gamma simulation. (F) The concentration of STAT1 significantly affected 
recruitment from STAT3 to (IFN-R-JAK)2*. 



independent of STATl disruption after IL-6 stimulation, 
although changing the STATS level significantly affected 
the formation of (IL6-gp80-gpl30-JAK)2*-STATl(data 
not shown). Therefore, the unbalanced competition be- 
tween STATl and STATS for IFNR and gplSO was not 
only the pivotal mechanism for the preferential activa- 
tion of IFN-gamma and IL-6, but it also determined the 
recruitment of STATl and STATS to the activated re- 
ceptor complexes. 

Responses of the crosstalk model to combined 
stimulation with IFN-gamma and IL-6 

We considered a combined stimulation protocol where 
the model was stimulated with IFN-gamma (0.1 nM) and 
IL-6 (0.1 nM) together for 12 h. Figure 4A shows that 



compared with the separate treatments, the combined 
stimulation induced higher activation of the JAK/STAT 
pathway. The dynamic responses of the JAK/STAT path- 
way were consistent with previous results reported in the 
literature [18]. In the previous section, we showed that 
IFN-gamma and IL-6 could activate both STATl and 
STATS, which may explain the higher activation of STATl 
and STATS after combined stimulation. However, when 
STATl could only be activated by IFN-gamma and STATS 
could only be activated by IL-6, combined stimulation still 
caused greater activation of STATl and STATS than sep- 
arate stimulation (Figure 4B). Thus, we inferred that other 
mechanism may play important roles in this phenomenon. 
PPl and PP2 are two different types of phosphatases that 
dephosphorylate STAT*s in the cytoplasm and nuclei. 
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Figure 4 Simulated time courses of species in the crosstalk model after combined stimulation with IFN-gamma and IL-6. (A) Combined 
stimulation witli IFN-gamma and iL-5 induced higlier activation of STAT1 and STATS than separate treatments. (B) A setting where STATl could 
only be activated by IFN-gamma and STAT3 could only be activated by IL-6 did not abolish the higher activation of STATl and STAT3 compared 
with individual stimulation experiment. (C) Combined stimulation with IFN-gamma and IL-6 reduced the level of PPl compared with individual 
treatments. (D) Combined stimulation with IFN-gamma and IL-6 reduced the level of PP2 compared with the individual treatments. (E) Without 
the formation of STATl/3 heterodimers, combined stimulation with IFN-gamma and IL-6 did not induce higher levels of (STATl N*)2 and 
(STAT3N*) compared with normal conditions. 



respectively [16,17]. In our model, the total amounts of 
both PPl and PP2 were fixed and only the monomers 
could combine and dephosphorylate the newly generated 
STAT*s. We then investigated the signal transduction pro- 
files of PPl and PP2 in response to IFN-gamma and/or 
IL-6 stimulation, and we found that combined stimulation 
with IFN-gamma and IL-6 could activate STATl and 
STATS, which resulted in lower levels of PPl and PP2 
than the separate treatments (Figure 4C-D). After com- 
bined stimulation, the activation of STATl prevented 
STATS from being dephosphorylated and vice versa. This 
mechanism contributed to the higher activation of STATl 
and STATS after combined stimulation with IFN-gamma 
and IL-6. 



Signal transduction via the JAK/STAT pathway 
depended on the formation of STAT homodimers [10], 
which are regarded as the main transcription factors 
during IFN-gamma and IL-6 signalling [2,8]. We further 
investigated whether combined stimulation with IFN- 
gamma and IL-6 could induce higher (STAT1N)''2 and 
(STAT3N)°'2 than separate treatments. Figure 4A shows 
that 12 h after combined stimulation, the STAT homodi- 
mers were not induced at a higher level than the separ- 
ate treatments. However, STATIN* -STATSN* reached 
their maximum concentration (200 nM) within about 
1 h, which was about S time higher than the individual 
treatment. The formation of STAT1N*-STATSN* greatly 
restricted the formation of STATs homodimers. After we 
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abolished the formation of STATl/3 heterodimers, the 
maximum concentrations of (STAT3N)*2 and (STAT1N)*2 
increased to about 100 nM with combined stimu- 
lation (Figure 4E). Combined stimulation with IFN- 
gamma and IL-6 led to greater activation of both STATl 
and STATS, but the formation of STATl/3 heterodimers 
played an important role in preventing mutual strengths 
between IFN-gamma and IL-6 signalling. 

Responses of the crosstalk model to successive IFN- 
gamma and IL-6 stimulation 

We analyzed previous studies that focused on the inter- 
actions between IFN-gamma and IL-6 signalling and 
found that their interactions were asymmetric. Bluyssen 
et al. reported that pre-treatment of EC with IFN- 
gamma significantly decreased STATS* induction by IL- 
6 without affecting the total amount of STATS [18]. By 
contrast, Kaur et al. reported that STATl activation 
induced by IFN-gamma was mainly unchanged after 
pre-treatment IL-6 or other gplSO-related cytokines in 
SH-SY5Y human neuroblastoma cells [24]. We tried to 
provide a reasonable explanation for the asymmetric 
interactions between IFN-gamma and IL-6 using simula- 
tion experiments with our model. 

First, we stimulated the model with IFN-gamma (0.1 
nM) for 12 h, which we started 2 h prior to IL-6 (0.1 nM) 
stimulation. IL-6 slighdy increases the level of STATS*, but 
pre-treatment with IFN-gamma significantly decreased 
STATS* induction by IL-6 (Figure 5A). This was consistent 
with the results reported by Bluyssen et al. [18]. SOCSS is 
a negative regulator of IL-6 signalling and it can be induced 
by IFN-gamma stimulation, so we deduced that SOCSS 
may have an important role during inhibition. When we 
knocked out SOCSS, the inhibitory effect of IFN-gamma 
on STATS* induction by IL-6 was eliminated completely 
(Figure 5B). These results showed that SOCSS was an es- 
sential component in the inhibition of IFN-gamma to IL-6 
signalling. We also found that the expression of SOCSS 
had a time delayed feedback, which significantiy increased 
1 h after IFN-gamma stimulation. Therefore, we deduced 
that temporal pre-treatment with IFN-gamma may not 
have induced sufficient SOCSS to inhibit IL-6 signalling. 
Figure 5C shows that temporal pre-treatment with IFN- 
gamma partiy inhibited IL-6 signalling and that the dur- 
ation of pre-treatment with IFN-gamma needed to be 
longer than 1 h to achieve this inhibition. We then investi- 
gated how pre-treatment with IL-6 affected the IFN- 
gamma signal response. Our simulation results showed 
that pre-treatment with IL-6 for 2 h only slightly reduced 
the amount of STATl* and did not inhibit the signal re- 
sponse of IFN-gamma (Figure 5D), while changing the dur- 
ation of the pre-treatment with IL-6 still had no obvious 
effect on the signal response of IFN-gamma. Moreover, 
pre-treatment for less than 1 h had almost no effect on the 



state of STATl* (Figure 5E). These simulation results were 
consistent with the results reported by Kaur et al. [24]. We 
inferred that the asymmetric interactions between IFN- 
gamma and IL-6 signalling were related mainly to the dif- 
ferent inhibition efficiencies of SOCSl and SOCSS. SOCSl 
could be induced by (STAT1N)*2 after IL-6 stimulation, 
but SOCSl induction by IL-6 is not sufficient to inhibit 
IFN-gamma signalling. After IFN-gamma stimulation, 
however, SOCSS could be induced by (STAT1N)*2 and 
(STATSN)*2, which achieved the inhibition from IFN- 
gamma to IL-6. 

The formation of STATl/S heterodimers also contribu- 
ted to the asymmetric interactions between IFN-gamma 
and IL-6 signalling. As previously described, the formation 
of STATl/S heterodimers enhanced the preferential signal 
transduction of IFN-gamma and IL-6 by sequestering a 
fraction of STATl* and STATS*. After we abolished the 
formation of STATl/S heterodimers, the maximum con- 
centration of SOCSl induced by IL-6 increased to 4.2 nM 
within about 1.5 h, and IL-6 exhibited a greater capacity 
for inhibiting IFN-gamma signalling (Figure 5F). Abolish- 
ing the formation of STATl/S heterodimers also enhanced 
the inhibition from IFN-gamma to IL-6 (Figure 5G). In 
addition, the mechanism of (STAT1N*)2 inducing SOCSS 
also played an important role in the asymmetric interac- 
tions. The concentration of (STATSN*)2 induced by IFN- 
gamma stimulation was very low due to the sequestering 
effect of STATl/S heterodimers. Therefore, we deduced 
that SOCSS induction by (STATSN*)2 was not sufficient 
to achieve the inhibition from IFN-gamma to IL-6. In- 
deed, when we abolished the (STAT1N)*2 induction of 
SOCSS, the inhibition from IFN-gamma to IL-6 was 
clearly mitigated (Figure 5H). IFN-gamma only slightly 
reduced the activation of STATS induction by IL-6, which 
did not agree with previous experimental observations 
[18]. Therefore, (STAT1N*)2 induction of SOCSS was es- 
sential for the inhibition from IFN-gamma to IL-6 during 
signal transduction. 

Argument for the non-competitive model 

Ho et al. reported that STATS did not affect the activation of 
STATl during type I interferon responses [2S]. They argued 
that STATS did not suppress STATl activation via tyrosine 
phosporylation and they excluded inhibitory mechanisms 
such as the competition for docking sites and the inhibition 
of signalling events upstream of STATl tyrosine phosphoryl- 
ation. Based on their argument, we proposed an alternative 
hypothesis for the combined pattern between STATs and the 
receptors for IFN-gamma and IL-6 in addition to our previ- 
ous competition model, which we refer to as the non- 
competitive model. In the non-competitive model, STATl 
and STATS do not compete for the same docking sites of 
IFNR and gplSO, so we considered that STATl and STATS 
could bind different phosphorylated docking sites in IFNR 
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Figure 5 Simulated time course of species in the crosstalk model of IL-6 after pre-treatment with IFN-gamma or IFN-gamma after pre- 
treatment with IL-6. (A) Pre-treatment with IFN-gamma significantly decreased the activation of STAT3 induced by iL-6. (B) Without S0CS3, the 
inhibitory effect of IFN-gamma on STAT3* induction by IL-5 was completely eliminated. (C) Temporal pre-treatment with IFN-gamma partly 
inhibited IL-6 signalling. (D) Pre-treatment with IL-6 reduced the amount of STATl* only slightly and did not inhibit the response of IFN-gamma. 
(E) Changing the length of pre-treatment with IL-6 did not affect the signal response of IFN-gamma. (F) Without the formation of STAT1/3 
heterodimers, IL-6 had a higher capacity to inhibit IFN-gamma signal transduction. (G) Without the formation of ST ATI/3 heterodimers, IFN- 
gamma had a higher capacity to inhibit IL-6 signal transduction. (H) Abolishing (STATl N)*2 induction of S0CS3 inhibited the IFN-gamma effects 
on IL-6. 
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and gpl30 via their SH2 domain. Based on these consid- 
erations, we added the new biochemical reactions (N17) - 
(N20) (see "non-competitive model" in Additional file 1) to 
simulate the mutually independent combination of STATl 
and STATS with the receptor complexes. STAT3 did not 
suppress the activation of STATl via tyrosine phosphoryl- 
ation, so we considered that the combination of STATS 
with receptor complexes did not suppress the phosphoryl- 
ation processes of STATl and vice versa. Based on these 
considerations, we added the new biochemical reactions 
(N21) - (N24) to simulate the mutually independent phos- 
phorylation of STATl and STATS within receptor com- 
plexes. In addition, we estimated the parameters for 
reactions (N17)-(N24) and re-estimated the parameters for 
reactions (N1)-(N6) to fit the non-competitive model 
(Additional file 1: Table S5). Detailed descriptions of the 
non-competitive model written in COPSI (using the file- 
name extension ".cps") are available in Additional file 2. 

First, we conducted simulation experiments using IFN- 
gamma and IL-6 stimulation separately, as in the previous 
study, and we aimed to estimate and optimize the parameters 
of the new non-competitive model. As shown in Figure 6A-B, 
our computational simulation indicated that the activation of 
STATl was much greater than STATS after IFN-gamma 
stimulation, whereas the activation of STATS was much 
greater than STATl after IL-6 stimulation, which agreed with 
previous experiment results despite some deviations in the sig- 
nal strength and durations [19,21]. Next, the same kinetic af- 
finities were used for IFNR and gplSO with STATl and 
STATS, respectively. IFN-gamma and IL-6 stimulation caused 
similar strong activation of STATl, STATS, SOCSl and 
SOCSS (Figure 6C-D), which did not agree with previous ex- 
perimental observations. These results demonstrated that the 
non-competitive model using our estimated parameter set 
could satisfy the basic simulation requirements and that it 
simulated the preferential activation of IFN-gamma and IL-6. 

Next, we executed the same series simulation experi- 
ments described in the previous sections using the non- 
competitive model. We found that the non-competitive 
model approximately simulated the crosstalk between IFN- 
gamma and IL-6 signalling in all conditions (data not 
shown), with the exception of the disrupted STATl and 
STATS conditions. Figure 6E shows that when we fixed the 
concentration of STATl at 1000 nM and changed the con- 
centration of STATS from 0, 500, 1000, to 2000 nM, the 
states of STATl* were obviously dependent on the level of 
STATS after IFN-gamma stimulation. In particular, STATl 
was more strongly phosphorylated, and for longer, with 
higher levels of STATS, which were completely different 
from the experimental observations reported by Regis et al. 
[22]. Similarly, when we fixed the concentration of STATS 
at 1000 nM and changed the concentration of STATl from 
0, 500, 1000, to 2000 nM, the states of STATS* were also 
obviously dependent on the level of STATl after IL-6 



stimulation (Figure 6F). Moreover, STATS was also more 
strongly phosphorylated and for longer with higher levels 
of STATl, which did not agree with the experimental 
observations of Dimberg et al. [57]. Indeed, the activation 
of STATs was determined by the swing between phosphor- 
ylation and dephosphorylation processes [6]. In the non- 
competitive model, STATl and STATS did not suppress 
the phosphorylation of each other, although they share the 
phosphatases PPl and PP2, so the activation of STATl pre- 
vented STATS from being dephosphorylated and vice 
versa. The disruption of the dynamic balance between the 
phosphorylation and dephosphorylation of STATl and 
STATS meant the responses of the non-competitive model 
were hypersensitive to the initial concentrations of STATl 
and STATS. Therefore, the non-competitive model may 
not reflect the crosstalk mechanisms between IFN-gamma 
and IL-6 signalling in a physiological environment. Thus, 
alternative hypotheses should be proposed based on new 
experimental observations. Here, we provide an integrated 
platform that facilitates the verification of possible hypoth- 
eses related to the crosstalk between IFN-gamma and IL-6 
signalling in future work. 

Sensitivity analysis of the competition model 

In this study, we used computational simulation to demon- 
strate the effectiveness of the competition model. Next, we 
determined the dynamic characteristics of our model by 
sensitivity analysis. The calculations used in the sensitivity 
analysis are shown in the "Methods" section. The upstream 
key regulatory molecules of (STAT1N*)2 and (STATSN*)2 
comprise cytokines, receptors, JAK, STATIC, STATSC, 
PPl, PP2, and SHP2, and their initial concentrations were 
considered in the sensitivity analysis. First, we applied a 
sensitivity analysis to determine the relative sensitivity of 
these components using IFN-gamma stimulation as the 
input. 

The results indicated that IFN-gamma, the receptor, JAK, 
STATIC, and PP2 were relatively important components of 
the IFN-gamma-induced JAK/STAT signalling pathway 
(Figure 7A). It was also noted that some of the initial con- 
centrations, such as those of PPl and STATSC, had littie 
impact on the time course of (STAT1N*)2. In a previous 
simulation, we found that the decline in the concentration 
of PP2 was more obvious than that of PPl during IFN- 
gamma and IL-6 signal transduction (Figure 4C-D). We 
deduced that the basic level of PPl exceeded the demand of 
IFN-gamma transduction, which agreed with the conclu- 
sions of the previous study. The concentration of STATS 
did not affect the activation of STATl, as described in the 
previous section. Moreover, we found that the parameters 
related to the key signalling components identified had rela- 
tively significant effects on the time course of (STAT1N*)2, 
as shown in Figure 7B. The reactions with the most sensitive 
kinetic parameters were STATIC phosphorylation by the 
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receptor complex of IFN-gamma and PP2 negative regula- 
tion. The SOCSl synthesis, degradation, and negative regu- 
lation processes were also very sensitive, which confirmed 
the important role of the negative feedback factor SOCSl 
during IFN-gamma signal transduction. 



Next, we performed a sensitivity analysis using IL-6 
stimulation as the input. The sensitivities of the initial con- 
centrations and the critical kinetic parameters are shown 
in Figure 7C and Figure 7D. The results indicated that the 
reactions with PP2 negative regulation and JAK 
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Figure 6 Simulated time courses of species in the non-competitive model with continuous exposure to IFNy (0.1 nM) or IL-6 (0.1 nM). 

(A) IFNy stimulation induced rnudi stronger activation of STATl tlian STATS. (B) IL-5 simulation induced much stronger activation of STAT3 than 
STATl. (C) With balanced kinetic affinities, IFNy stimulation induced similar activation of STATl and STATS. (D) With balanced kinetic affinities, IL-6 
stimulation induced similar activation of STATl and STATS. (E) STATl was phosphorylated more strongly and for longer with a higher level of 
STATS after IFNy simulation. (F) STATS was phosphorylated more strongly and for longer with a higher level of STATl after IL-6 simulation. 
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Figure 7 Sensitivity analysis of the competition model. (A, B) indicate tlie sensitivities of the initial molecule concentrations and the reaction 
rate constants after IFN-gamma stimulation. (C, D) indicate the sensitivities of the initial molecule concentrations and the reaction rate constants 
after IL-5 simulation. The definition of the parameter sensitivity can refer to the "Methods" section. The positive (or negative) sensitivity values 
indicate consistent (or opposite) percentage changes in output to the perturbation in the parameter compared to the nominal solution. The red 
and blue bars represent positive and negative sensitivity values, respectively. For brevity, we only show the top 20 parameters with high 
sensitivity values. 



combination were highly sensitive to IL-6 stimulation. As 
described in previous sections, the concentration of STATl 
did not significantly affect the activation of STATS after 
IL-6 stimulation, which was also confirmed by the sensitiv- 
ity analysis. Similarly, we found that the SOCS3 synthesis, 
degradation, and negative regulation processes were highly 
sensitive to IL-6 stimulation, which confirmed the import- 
ant role of the negative feedback factor SOCS3 during IL- 
6signal transduction. In particular, we found that the phos- 
phorylation of STATS in STATl/S heterdimers in the nu- 
cleus also had a high impact on the state of (STAT3N*)2. 
Overall, the sensitivity analysis determined the sensitive 
components and parameters during JAK/STAT signal 
transduction. These results provide valuable information 
that may inform further investigations of regulation and 
drug target identification in aberrant pathways. The 
detailed results of the sensitivity analysis of the competition 
model are shown in Additional file 1: Tables S6-S7. 

Sensitivity analysis of the non-competitive model 

We also investigated the dynamic characteristic of the 
non-competitive model by sensitivity analysis. First, we 
applied sensitivity analysis to determine the critical com- 
ponents with dominant effects where IFN-gamma stimu- 
lation was used as the input. Figure 8A shows the 
sensitivities to the initial concentrations, while the sensi- 
tivities of the critical kinetic parameters are shown in 
Figure 8B. There were some deviations in magnitude but 



the non-competitive model yielded similar sensitivity results 
to the competition model after IFN-gamma stimulation, ex- 
cept the concentration of STATSC had a higher sensitivity 
to (STAT1N*)2 in the non-competitive model. This con- 
firmed our previous conclusion that the responses of IFN- 
gamma in the non-competitive model were hypersensitive 
to the initial concentrations of STATS. Next, we performed 
a sensitivity analysis using IL-6 stimulation as the input. 
Figure 8C shows the sensitivities to the initial concentra- 
tions, while the sensitivities to the critical kinetic parameters 
are shown in Figure 8D. We also found that the concentra- 
tion of STATIC was highly sensitive to (STAT3N*)2 in the 
non-competitive model. In the previous section, simulation 
results demonstrated the non-competitive model could not 
accurately reflect the signal transduction of IL-6 and IFN- 
gamma under the condition of disrupting of STATs. Here, 
the sensitivity analysis confirmed the high sensitivity of 
STATs in the non-competitive model. The detailed results 
of the sensitivity analysis of the non-competitive model are 
shown in Additional file 1: Tables S8-S9. 

Discussion 

The results were affected by the limitations of our model. It 
is important to consider these limitations because they form 
the basis for fiiture improvements. First, we only considered 
three possible levels of crosstalk between IFNy and IL-6 sys- 
tems. Haan et al. showed that, after IL-6 stimulation, STATl 
phosphorylation was entirely dependent on JAKl kinase 
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Figure 8 Sensitivity analysis of the non-competitive model. (A, B) indicate the sensitivities of the initial nnolecule concentrations and the 
reaction rate constants after IFN-gamma stimulation. (C, D) indicate the sensitivities of the initial molecule concentrations and the reaction rate 
constants after IL-6 simulation, respectively. The definition of the parameter sensitivity can refer to the "Methods" section. The positive (or 
negative) sensitivity values indicate consistent (or opposite) percentage changes in output to the perturbation in the parameter compared to the 
nominal solution. The red and blue bars represent positive and negative sensitivity values, respectively. For brevity, we only show the top 20 
parameters with high sensitivity values. 



activity whereas STATS activation was not [48]. Qing et al. 
showed that, in response to IFN-gamma, SCR-famUy kinases 
had to activate STATS (but not STATl) via tyrosine phos- 
phorylation [19]. Some mechanisms that contribute to the 
specific signal responses of IFN-gamma and IL-6 were not 
considered in this study. Second, our computational model 
was a simplification of a biological process and it might not 
reflect the true regulatory mechanism for the signal interac- 
tions. Third, the simulation results in this study were affected 
by assumptions made based on our experience, e.g., the 
STAT 1/3 heterodimers could be dephosphorylated by PPl 
and PP2 in the cytoplasm and nucleus, respectively. Fourth, 
our model had very rich dynamics parameters that were not 
fiJly explored in this study. For simplicity, protein turnover, 
receptor recycling, de novo synthesis and the degradation of 
transcription factors were not included in the current model. 
The ultimate aim of our research is to construct a universal 
model that accurately reflects the crosstalk between IFN- 
gamma and IL-6 signals. However, the signal responses to 
cytokines may be cell type-dependent. Caldenhoven et al. 
reported that the activation of STATS induced by IFN- 
gamma was lineage-specific in human neutrophils [58]. 
Zhang et al. showed that IL-6 stimulation could induce 
STATl phosphorylation in a dose- and time-dependent 
manner in Ml, R2 and U937 cells, although it had litde effect 
on STATl phosphorylation in 7TD1 and TFl cells [21]. 
Bluyssen et al. argued that IL-6 did not activate STATl in 
EC [18]. Furthermore, the concentrations of molecules 



within the JAK/STAT pathway, such as STATs, are cell type- 
dependent [6]. The experimental data are limited and not 
systematic, so we had to create our model based on experi- 
mental observations of many diSferent cell types, while we 
neglected certain contradictory experimental observations. 
Therefore, it is not possible to expect that the dynamics pre- 
dicted by our model will apply universally to all types of cells. 
Thus, structure and parameters of the model may need some 
adjustment to reflect signal transduction by IFN-gamma and 
IL-6 in certain cell types. Finally, our crosstalk model was 
based on experimentally established interactions, but further 
experimental verification and improvements are required. 

Our simulation results showed that STATl/S heterodi- 
mers have three important functions during signal trans- 
ductions from IFN-gamma to IL-6. First, the formation 
of STATl/S heterodimers enhances the preferential sig- 
nal transduction by IFN-gamma and IL-6 because it 
sequesters a fraction of STATl* and STATS*. Second, the 
formation of STATl/3 heterodimers prevents mutual 
reinforcement between IFN-gamma and IL-6 signalling. 
Finally, the formation of STATl/S heterodimers limits the 
reciprocal inhibition of IFN-gamma and IL-6 signalling. In 
our simulations, therefore, the formation of STATl/3 het- 
erodimers dramatically affected the interaction between 
the IFN-gamma and IL-6 systems, which suggests that 
STATl/S heterodimers may be a potential target for recti- 
fying abnormal signal transduction by IFN-gamma and 
IL-6. The functional interference of STATS homodimers 
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using STATS transcription decoys or small molecules in 
structure-activity relationship (SAR) studies could suc- 
cessfully inhibit the growth of tumour cells [59,60]. How- 
ever, the therapeutic potential of altering the formation of 
STATl/3 heterodimers has not been fully investigated. 
Thus, further research is still required. 

This is the first effort to construct a mathematical 
model of the crosstalk between IFN and IL-6 signal- 
ling. Moreover, our simulation results and theoretical 
findings provide new insights into the dynamical 
integration of IFN and IL-6 signals. The lack of ex- 
perimental data and our current superficial under- 
standing of signal transduction mean there are still 
many defects in our crosstalk model. We will follow 
current studies related to this issue and improve the 
quality and veracity of model in our future research. 

Conclusions 

Based on previous models and new experimental obser- 
vations, we developed the first crosstalk model of IFN- 
gamma and IL-6 signalling. This theoretical study suc- 
cessfully reproduced key experimental findings and 
reached some definitive conclusions. First, the unba- 
lanced competition between STATl and STATS for 
IFNR and gplSO led to preferential activation of IFN- 
gamma and IL-6. At the same time, the formation of 
STATl/S heterodimers enhanced preferential signal 
transduction by sequestering a fraction of STATl* and 
STATS*. Moreover, SOCSs with SHP-2 limited the con- 
centration of the activated receptor complexes of IFN- 
gamma and IL-6, which also contributed to the preferen- 
tial activation of IFN-gamma and IL-6. Second, the 
unbalanced competition between STATS and STATl 
was the pivotal mechanism during the mutual switch be- 
tween IFN-gamma and IL-6 signals after knocking out 
STATl or STATS. Finally, the formation of STATl/S 
heterodimers prevented the mutual reinforcement of 
IFN-gamma and IL-6 signalling, and also limited the re- 
ciprocal inhibition between IFN-gamma and IL-6 signal- 
ling. Moreover, the process of STATl homodimer 
induction of SOCSS also contributed to asymmetric 
interactions between IFN-gamma and IL-6 signalling. 

Methods 

Three classes of kinetics were involved in our crosstalk 
model. The mass-action equation for the molecular 
combination and decomposition was as follows: 

A + Bi^C + D 

where /<f and are the rates of the forward and reverse 
reactions. The Michaelis-Menten function (actually be 



split into two Mass-Action reactions) for the enzymatic 
reactions was as follows: 

E + S:f±ES^E + P 

where E is the enzyme; S is the substrate; £5 is the 
enzyme-substrate complex; P is the product; and Ki, /<2 
and /<3 are the rates of the reactions. The equation for 
substances undergoing translocation was as follows: 



where Ac represents a molecule in the cytoplasm; Aj^ is 
Ac translocated from the cytoplasm to nuclei; B is the 
catalyst; and Kj is the rate of the reaction. 

Based on Yamada et al. [28] and Moya et al. [SI], we 
established the components of the IFN-gamma and IL-6 
pathways, respectively, where the model parameters of 
the two components were left unchanged. The equations 
of these reactions are shown in Additional file 1 and the 
kinetic parameters of these reactions are provided in 
Additional file 1: Table S2-SS. Based on the possible 
physiological mechanisms of IFN-gamma and IL-6 
cross-regulation, we also added 16 new reactions (Nl)- 
(N16) in our crosstalk model. The equations of these 
reactions are listed below (New biochemical reactions 
added to the crosstalk model) while a graphic descrip- 
tion of these reaction is shown in Figure 9. The kine- 
tic parameters of these new reactions are provided in 
Additional file 1: Table SI. 

New biochemical reactions added to the crosstalk model 

STATSC binding: this positive reaction describes 
STATSC binding with the receptor complex of IFN- 
gamma while the reverse reaction describes the dissoci- 
ation of STATSC. 

{IFN-R-/AK)2 

* +STAT3C% {IFN -R- JAK)2 

* -STAT?,C 

(ATI) 

STATSC activation: this reaction describes the activa- 
tion of STATSC by the receptor complex of IFN-gamma. 

{IFN -R-JAK)2 

* - STA T3C ^ {IFN -R- JAK)2 

* + STAT?,C* 

{N2) 

STATSC* binding: this positive reaction describes 
STATSC* binding with the receptor complex of IFN- 
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Figure 9 Schematic description of the newly added 
biochemical reactions in the crosstalk model. In the diagram: R 
represents a receptor; C represents a species in the cytoplasm; N 
represents a species in the nucleus; 2 represents dimers; and * 
represents the activation state. mRNA-S0CS3N represents the S0CS3 
mRNA in the nucleus. Arrows indicate the reaction direction. 



gamma while the reverse reaction describes the dissoci- 
ation of STAT3C*. 



{IFN -R- JAK)2 * ^STAT-iC 

* {IFN -R- JAK)2 * -STAT?,C* 



(N3) 



STATIC binding: this positive reaction describes 
STATIC binding with the receptor complex of IL-6 
while the reverse reaction describes the dissociation of 
STATIC. 

(/L6 - gpS,0 - gpUQ - JAK)2 

* +STATIC ?=± (716 - gp80 - gpUQ - JAK)2 

* -STATIC 

(N4) 

STATIC activation: this reaction describes the activa- 
tion of STAT3C by the receptor complex of IL-6. 

{IL6 - gp80 - gpl30 - JAK)2 

* -static"-^ {IL6 - gp80 - gpl30 - JAK)2 

* +STATIC* 

(N5) 

STATIC* binding: this positive reaction describes 
STAT3C* binding with the receptor complex of IL-6 



while the reverse reaction describes the dissociation of 
STATIC*. 



{IL6 - gpS,0 - gpno - JAK)2 * +STATIC 

* {IL6 - gp80 - gpUQ - JAK)2 

Kfc-6 

* -STATIC* 



(N6) 



mRNA-SOCS3 transcript: this reaction describes 
SOCS3 mRNA transcription after its induction by 
STATl homodimers in the nucleus. 

d[mRNA - S0CS3N]/dt 

= V max{STATlN*)2/[{STATlN*)2 + Km] 

(N7) 

STATl/3 heterodimer formation: this positive reaction 
describes the formation of STATl/3 heterodimers in the 
cytoplasm while the reverse reaction describes their dis- 
sociation. 

STATIC * +STAT3C * <^ STATIC 

* -STAT3C* (N8) 

STATl/3 heterodimer formation: this positive reaction 
describes the formation of STATl/3 heterodimers in the 
nucleus while the reverse reaction describes their dis- 
sociation. 



«if9 

STATIN * +STAT3N * STATIN 

ntr9 

* -STAT3N'i^ 



(N9) 



STATl/3 heterodimer translocation: this reaction describes 
the nuclear importation of STATl/3 heterodimers. 



STATIC * -STAT3C * 
* -STAT3N* 



STATIN 



(NIO) 



PPl binding: this positive reaction describes PPl bind- 
ing with STATl/3 heterodimers in the cytoplasm while 
the reverse reaction describes their dissociation. 

PPl + STATIC * -STAT3C * ^ PPl 

"fell 

-STATIC * -STAT3C* 

(Nil) 

PP2 binding: this positive reaction describes PP2 bind- 
ing with STATl/3 heterodimers in the nucleus while the 
reverse reaction describes their dissociation. 

PP2 + STATIN * -STAT3N * ^ PP2 

Ufa. 12 

-STATIN * -STAT3N* 

(N12) 
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STATIC inactivation: this reaction describes the in- 
activation of STATl in STATl/3 heterodimers by PPl in 
the cytoplasm. 

PPl - STATIC * -STAT3C * PPl 

+ STATIC + STAT'iC'i' (N13) 

STAT3C inactivation: this reaction describes the inactivation 
of STATS in STATl/3 heterodimers by PPl in the cytoplasm. 

PPl - STATIC * -STAT3C * PPl 

+ STATIC * +STAT3C (N14) 

STATIN inactivation: this reaction describes the in- 
activation of STATl in STATl/3 heterodimers by PP2 in 
the nucleus. 

PPl - STATIN * -STAT3N * PPl 

+ STATIN + STAT3N* (N15) 

STAT3N inactivation: this reaction describes the inactivation 
of STATS in STATl/S heterodimers by PP2 in the nucleus. 

PPl - STATIN * -STAT3N * PPl 

+ STATIN * +STAT3N (N16) 

where R represents the receptor; C represents a species within 
the cytoplasm; N represents a species within nuclei; 2 repre- 
sents dimers; and * represents the activation state. mRNA- 
SOCS3N represents the SOCSS mRNA in the nucleus; 
(STAT1N*)2 represents phosphorylated STATl homodimers 
in the nucleus; V^ax represents the maximum reaction 
velocity of the SOCS3-mRNA transcription; represents 
half of the substrate concentration when the enzymatic reac- 
tion reaches its maximum reaction velocity. 

Using this reaction scenario, we established systematic 
models that were described using ordinary differential 
equations (ODEs). COPASI (http://www.copasi.org/) and 
CellDesigner (http://www.celldesigner.org/) were used to 
calculate the model states and solve the ODE models 
[61,62]. Detailed descriptions of the model written in 
COPSI (using the filename extension ".cps") are available 
in Additional file 2. 

Sensitivity analysis 

The sensitivity analysis determined the change in the 
output after a tiny perturbation in the initial state or the 
reaction parameters. Thus, a perturbation Ax in the par- 
ameter set induced an overall state change Ay and the 
sensitivity was defined as Ay I Ax, which was normalized 
with respect to the actual values of the phenomenon 
investigated and the parameter x/y. Experimental studies 
have shown that phosphorylated STAT dimers in the 
nucleus are critical transcription factors for target gene 
activation in the JAK/STAT pathway. Therefore, our 



sensitivity analysis considered (STAT1N*)2 and 
(STATSN*)2 as the outputs of IFN-gamma and IL-6 sig- 
nalling, respectively. In this study. Ay was defined as the 
change in the concentration of (STAT1N*)2 or 
(STATSN*)2 after 1.5 h, which adequately represented the 
signal transduction strength of IFN-gamma or IL-6. The 
perturbation ranges of each parameter were set as 30% of 
the actual absolute value to avoid errors in the sensitivity 
estimates. COPASI and CellDesigner were used to calcu- 
late the sensitivity. 
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Additional file 1: This file contains one figure, nine tables, and 
description of the biochemical reactions in the crosstalk model and 
the non-competitive model. Figure SI shows the overall structure of 
the crosstalk model. Table SI lists the parameters and values of the 
newly added reactions in the crosstalk model. Table S2 lists the 
parameters and values of the IFN-gamma part of the crosstalk model. 
Table S3 lists the parameters and values of the IL-6 part of the crosstalk 
model. Table S4 lists the state variables of the model and their initial 
values. Table S5 lists the re-estimated parameters and their values in the 
non-competitive model. Table S6 lists the results of the sensitivity 
analysis for variations in the concentrations of the pathway components 
in the competition model. Table S7 lists the results of the sensitivity 
analysis for variations in the kinetic parameters in the competition model. 
Table S8 lists the results of the sensitivity analysis for variations in the 
concentrations of pathway components in the non-competitive model. 
Table S9 lists the results of the sensitivity analysis for variations in the 
kinetic parameters in the non-competitive model. All of the biochemical 
reactions investigated in this study are described in this file, which 
contains the newly added biochemical reactions, the IFN-gamma and IL- 
6 parts of the crosstalk model, and the newly added reactions in the 
non-competitive model. 

Additional file 2: This file contains the crosstalk model and the 
non-competitive model written in COPASI, which can be executed 
using COPASI version 4.7. 
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